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Abstract 

Despite the importance of dimetliylsulplionlopropionate (DMSP) in tine global sulphur cycle and climate regulation, the 
biological pathways underpinning its synthesis in marine phytoplankton remain poorly understood. The intracellular 
concentration of DMSP increases with increased salinity, increased light intensity and nitrogen starvation in the diatom 
Thalassiosira pseudonana. We used these conditions to investigate DMSP synthesis at the cellular level via analysis of 
enzyme activity, gene expression and proteome comparison. The activity of the key sulphur assimilatory enzyme, adenosine 
5'-phosphosulphate reductase was not coordinated with increasing intracellular DMSP concentration. Under all three 
treatments coordination in the expression of sulphur assimilation genes was limited to increases in sulphite reductase 
transcripts. Similarly, proteomic 2D gel analysis only revealed an increase in phosphoenolpyruvate carboxylase following 
increases in DMSP concentration. Our findings suggest that increased sulphur assimilation might not be required for 
increased DMSP synthesis, instead the availability of carbon and nitrogen substrates may be important in the regulation of 
this pathway. This contrasts with the regulation of sulphur metabolism in higher plants, which generally involves up- 
regulation of several sulphur assimilatory enzymes. In 7. pseudonana changes relating to sulphur metabolism were specific 
to the individual treatments and, given that little coordination was seen in transcript and protein responses across the three 
growth conditions, different patterns of regulation might be responsible for the increase in DMSP concentration seen under 
each treatment. 
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Introduction 

Marine phytoplankton play a key role in the global sulphur 
cycle through the synthesis of dimethylsulphoniopropionate 
(DMSP), the major precursor of the volatile sulphur compound 
dimethylsulphide (DMS). DMS transfers sulphur from the oceans, 
which are a major sulphur reservoir, to the relatively sulphur- 
limited land [1]. Furthermore, once in the atmosphere, DMS 
oxidises to form aerosol particles and thereby contributes to 
cooling the climate, directly through the reflection of solar 
radiation and indirectiy through the formation of cloud conden- 
sation nuclei [2] . The global annual flux of DMS from the oceans 
into the marine atmosphere is estimated at be between 15 and 
33 Tg sulphur per year [3] and Gunson et al. [4] used a modelling 
approach and scenario to demonstrate that halving the DMS flux 
could have a significant effect on radiative forcing, increasing 
surface temperatures by around 1.6°C. More recent debate 
highhghts additional phytoplankton-derived compounds that 



may also produce particles and reflect radiation from the Sun 
back into space [5] . 

Considerable variability in DMSP production has been 
observed between phytoplankton taxa [6]. The Dinophyceae 
and the Prymnesiophyceae are the highest DMSP producers with 
intracellular concentrations in excess of several hundred mmol 
1—1 in some species, although diatoms and members of other 
groups can also produce significant amounts. The chlorophytes on 
the other hand produce very littie DMSP [6], with the exception of 
seaweeds such as Ulva lactuca [7]. In addition to taxonomy, 
external factors, such as nutrient availability, salinity and 
temperature can cause variations in intraceUular DMSP concen- 
tration. Multiple cellular roles have been proposed for DMSP, 
including as an osmolyte [7,8], cryoprotectant [9], grazing 
deterrent [10], antioxidant defence [11] and as an 'overflow 
metabolite' in dissipation of excess energy, carbon and reducing 
equivalents when growth and photosynthesis are unbalanced [12]. 
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The primary role of DMSP in the cell, however, remains uncertain 
and might vary within the lifespan of a cell and between species. 

DMSP is generated from methionine [13] and Gage et al. [14] 
used ''^S-labelling to determine the steps converting this amino 
acid to DMSP in the green macroalga Ulva intestinalis. Methionine 
first undergoes transamination to 4-methylthio-2-oxobutyrate 
(MTOB) and then reduction to 4-methylthio-2-hydroxybutyrate 
(MTHB), before 5'-methylation to 4-dimethylsulphonio-2-hydroxy- 
butyrate (DMSHB) and finally oxidative decarboxylation to 
DMSP (Figure 1) [14]. The key intermediate, DMSHB, has also 
been identified in tlu- marine microalgae Emiliania huxleyi, 
Tetaselmis sp. and the diatom Melosira nummuloides, indicating that 
the same pathway exists in diverse algal classes [14]. The activities 
of putative, substrate-specific enzymes for the first three steps of 
this pathway were measured in U. intestinalis [15], although the 
enzymes have not been identified. 

Despite the importance of DMSP production in the sulphur 
cycle and climate regulation, and the availability of genome 
sequencers for s(;veral DMSP producing phytoplankton, there has 
been limited research into basic sulphur metabolism in this group. 
We have very little direct knowledge of this pathway in most 
phytoplankton species and know almost nothing about its 
regulation. A detailed understanding of the biological mechanism 
of DMSP production is required to improve our ability to predict 



Sulphur 
assimilation 



Serine 
Serine acetyl- I 
transferase 
(SAT) I 
0-acetylserine 



Sulphate (external) 

Sulptiate transporter ^ 

Sulphate (internal) 

ATP sulphurylase ± 
(ATPS) ^ 
APS 

APS reductase 

(APR) r 

Sulphite 
Sulphite reductase 
(SIR) 

Sulphide 
OAS thiol lyase\ 
JOASTL) 



'I 



1 

_Cireteine_ 



Aspartate 
I 

O-phosphohomoserine- 

; 

Threonine 



Methionine 
synthesis 



Cystathionine y-synthase 
(CgS) 



-► Cystathionine 

Cystathionine p-lyase 
▼ (CbL) 
Homocysteine 



1 Methionine synthase 
Methionine 



TRANSAMINATION 

4-methylthio-2-oxobutyrate 

\ REDUCTION 
4-methylthio-2-hydroxybutyrate 

^ S-METHYLATION 

4-d i methylthio-2-hydroxybutyrate 

\ OXIDATIVE 
DMSP DECARBOXYLATION 



synthesis 



Figure 1. A proposed pathway of sulphate assimilation and 
DMSP biosynthesis in algae. Enzymes between the methionine and 
DIVISP have not yet been identified. The putative reactions proposed by 
[13] are shown in capitals. 
doi:1 0.1 371/journal.pone.0094795.g001 



how environmental factors might affect this process. Comparing 
cell protein and transcript abundance has the potential to offer 
new insight into the pathway of sulphur assimilation and DMSP 
biosynthesis, and its regulation, in marine phytoplankton. 

Diatoms are not generally considered major DMSP producers, 
but the synthesis of DMSP liy members of this group has been 
shown to be up-regulated under certain growth conditions [6,16]. 
For example, Bucciarelli and Sunda [16] measured increased 
DMSP concentration in the diatom Thalassiosira pseudonana under 
nitrat(', phosphate and silicate limitation. Given that diatoms are 
one of the most abundant groups of phytoplankton, accounting for 
approximately 20% global primary production, they might 
therefore make a greater contribution to DMSP production than 
suggested by early studies on actively growing, nutrient-replete 
batch cultures [6]. In addition, the abihty to up-regulate DMSP 
production is a useful characteristic for investigating the regulation 
of its synthesis pathway. T. pseudonana is also an excellent model for 
diatom biology because its genome has been sequenced [17] and 
there are various molecular tools for this species. Our aim was to 
use T. pseudonana to address the control of DMSP synthesis at a 
cellular level. 

We analysed the proteome response of T. pseudonana to 
increased salinity and increased light intensity, two conditions 
that we have confirmed to induce DMSP synthesis in this species, 
using 2-dimensional gel electrophoresis. We then compared these 
findings to our previous study on the proteome response of T. 
pseudonana to nitrogen starvation [18], a growth condition that also 
increases intrac:ellular DMSP concentration in this species [16]. 
The overall hypothesis was that enzymes related to DMSP 
synthesis would be found amongst proteins that increased or 
decreased in abundance across all three of these growth 
conditions. In addition, since DMSP is an important sink for 
reduced sulphur, we also compared the transcript abundance of 
genes involved in the uptake and assimilation of sulphate and 
measured the activity of the key enzyme of the pathway, adenosine 
5'-phosphosulphate reductase (APR) under increased salinity, 
increased light intensity and nitrogen starvation. 

Materials and Methods 

Culturing 

Axenic cultures of T. pseudonana (CCMP 1335, National Centre 
for Marine Algae and Microbiota, Bigelow, USA) were grown in 
batch culture in ESAW (enriched seawater, artificial water) 
medium [19] at 15oC with a 14:10 light:dark cycle. Unless 
otherwise stated, cultures were grown at 115 (tmol photons 
m— 2 s— 1 based on an immersed measurement with a Scalar 
PAR Irradiance Sensor Q_SL 2101 (Biospherical Instruments Inc., 
San Diego, USA). Cultures were regularly checked for bacterial 
contamination by 4',6-diamidino-2-phenyhndoledihydrochloride 
(DAPI) staining [20]. Cell number and volume were measured 
with a Beckman Coulter Multisizer 3 Analyser (Beckman Coulter 
Ltd, High Wycombe, UK) and variable to maximum fluorescence 
ratio (Fv/FJ with a Walz Phyto-Pam phytoplankton analyser 
(Heinz Walz GmbH, Effeltrich, Germany). Experiments were 
conducted in tripKcate for control and treated cultures. AH 
measurements and samples for further analysis were collected 
3 hours into the light period. 

Alteration of Salinity 

The salinity of ESAW growth medium was adapted by either 
increasing or decreasing the concentrations of all salts in the 
ESAW artificial seawater base recipe [19] with the exception of 
bicarbonate, because changes in its concentration affected culture 
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growth rate and yield. T. pseudonam cells, acclimated to a salinity of 
10 psu for a minimum of four rounds of subculture (approximately 
4 weeks) prior to the experiment, were grown in triplicate to ca. 
1x106 cells ml— 1 . At this point the cultures were divided equally 
and an equal volume of medium with a salinity of either 10 psu to 
maintain the salinity, or 60 psu, to achieve a final salinity of 35, 
was added to each portion. Samples were taken for transcript and 
proteome comparison 48 h after the salinity change, on day 4. 

The salinity was adjusted by the dilution of all salts because it is 
a more environmentally relevant treatment than only changing 
levels of sodium chloride, however, this also decreased the 
concentration of sulphate from 25 mM, as in standard ESAW 
medium, to 5 mM. This decrease in sulphate concentration 
without changing overall salinity was tested and found to have no 
effect on growth rate or intracellular DMSP concentration (Figure 
SI). 

Alteration of Light Intensity 

T. p.seudonana cultures acclimated to a reduced light intensity of 
50 |J,mol photons m— 2 s— 1 (achieved using a neutral density 
filter) for a minimum of four rounds of subculturing (approxi- 
mately 4 weeks) were grown to a density of ca. 5x1 05 cells ml— 1 . 
Cultures were then either kept at 50 (Xmol photons m— 2 s— 1 or 
exposed to a high light intensity of 1000 |J,mol photons m— 2 s— 1. 
Samples were taken for transcript and proteome comparison 48 h 
after cultures were exposed to increased light intensity, on day 8. 

Nitrogen Starvation 

As described in Hockin et al. [18] T. p.seudonana cultures were 
started with initial concentrations of 550 nM (standard ESAW) or 
30 |iM nitrate. Cultur(;s with an initial concentration of 30 |iM or 
550 nM nitrate became yield-limited at ca. 1 xl06 cell ml— land 
ca. 2x106 cell ml— 1, respectively. Samples were taken for 
transcript and proteome comparison on day 3 of the experiment. 

Intracellular DMSP Concentration 

Intracellular DMSP concentration was measured using head- 
space gas chromatography, developed by [21]. Depending on the 
cell density, 3 to 7 ml of culture was collected by filtration onto 
25 mm nucleopore track-etch membranes with a 1 |lm pore 
(Whatman) at a vacuum pressure no greater than 1 5 kPa. Samples 
were then handled as described in [22], with the exception that 
150 nl of headspace gas was sampled and injected into the GC. 
Technical triplicates were conducted for this method. 

Gene Expression Analysis 

T. pseudonam cells were collected by quick filtration of 250 ml 
culture aliquots onto 47 mm diameter 1 ^m nucleopore mem- 
branes (Whatman) at a vacuum pressure of 35 kPa. Filters were 
immediately frozen in liquid nitrogen and stored at — 80°C. Cells 
were washed from the membrane with 500 |J.l buffer (80 mM Tris 
pH 9, 5% SDS, 150 mM LiCl, 50 mM EDTA), and homoge- 
nised. Total RNA was extracted by phenol/chlorophorm and LiCl 
precipitation. 

Quantitative real-time RT-PCR (qPCR) was performed using 
gene-specific primers with efficiencies between 87 and 105% 
(Table SI) and the fluorescent intercalating dye SYBR Green 
(Sigma-Aldrich), in a DNA engine OPTICON2 continuous 
fluorescence detector (Bio-Rad). The programme steps were: 
2 min at 95°C, then 40 cycles consisting of 15 s at 95°C, 15 s at 
60°C, 30 s at 72°C and 10 min at 72°C, followed by a standard 
dissociation protocol to ensure that each amplicon was a single 
product. The relative difference in transcript abundance between 



the control and treatment samples was calculated using the AACt 
method [23]. Beta- tubulin (TUB 3) was used as the reference gene, 
and its stability was confirmed using geNorm software [24]. The 
RT-PCR reactions were performed in triplicate for each of the 
three independent biological replicates. 

APR Activity 

Depending on the cell density, 5 to 15 ml of culture was 
centrifuged and resuspended in 500 [i\ of extraction buffer 
(45.45 mM Na/KP04, 27 mM NazSOj, 0.45 mM AMP, and 
9. 1 mM DTE). The cells were disrupted by sonication (three times 
for 10 sec at ca. 15 microns on ice with a Soniprep 150 probe; 
MSE, London, UK) and APR activity was measured as the 
production of [35 Sjsulphite assayed as acid volatile radioactivity 
formed in the presence of [35 S] APS and dithioerythritol (DTE) as 
reductant [25]. APR activity was normalized to protein concen- 
tration, which was measured according to the Bradford assay [26] 
using a Bio-Rad Protein kit (Figure S2). 

Proteome comparison by 2-dimensional gel 
electrophoresis 

The cultures were collected by filtration of multiple 250 ml 
eiliquots onto, 47 mm diameter nucleopore membranes with a 
1 |im pore size (Whatman) at a vacuum pressure of 35 kPa. Filters 
were immediately frozen in liquid nitrogen and stored at — 80°C. 
For each of the three independent cultures, proteins from two 
filters were extracted, and 100 ig protein, as quantified with an 
Ettan 2-D Quant kit (GE Healthcare, Chalfont, UK) was 
subjected to 2-dimensional (2-D) gel electrophoresis and imaged 
exactly as described in [18]. Gel images were compared using 
Progenesis SameSpots analysis software (v4.1; Nonlinear Dynam- 
ics Ltd, Newcastie Upon Tyne, UK), which includes automatic 
background subtraction and normalization. Protein spots with 
altered levels of expression (1.5 fold change and q<0.05, t-test 
corrected for false discovery rate (FDR) using the Benjamini 
Hochberg procedure [27]) under treated versus control conditions, 
in all three growth conditions, were excised from the gel using a 
ProPick excision robot (Genomic Solutions). The excised proteins 
were then manually in-gel trypsin digested and analysed by 
peptide mass fingerprinting [28], as described in [18]. 

Results and Discussion 

Growth Conditions Affecting DMSP 

The identification of conditions that increase the concentration 
of DMSP in T. pseudomna cells was a critical first step in enabling 
us to study ceUs with up-regulated DMSP synthesis. We confirmed 
that the intraceUular DMSP concentration of T. pseudomna 
increases with increased salinity, increased light intensity and 
nitrogen starvation. CeUs transferred from 10 psu to 35 psu had 
an intracellular DMSP concentration of 8.9 mM after 48 hours 
whilst its concentration in ceUs maintained at 10 psu remained 
below the level of detection (Figure 2A,B). When samples were 
collected for proteome and transcriptional analysis, 48 h after the 
increase in salinity, no difference in Fv/Fm between the two 
treatment groups was observed (Figure 2C). Exposure to a high 
light intensit}' of 1000 )a.mol photons m— 2 s— 1 for 48 hours lead 
to a 8-fold higher intracellular DMSP concentration compared to 
cells maintained at 50 |Xmol photons m— 2 s— 1 (Figure 2D,E). 
The Fv/Fm of cultures exposed to a higher light intensity 
decreased within 24 h and remained below that of cultures 
maintained at a lower light intensits' at 48 h when samples were 
collected for proteome and transcriptional analysis (Figure 2F). 
The intraceUular DMSP concentration of nitrogen starved T. 



PLOS ONE I www.plosone.org 



3 



April 2014 | Volume 9 | Issue 4 | e94795 



Regulation of DMSP Synthesis 



pseudonana on day 3 was 3.5-fold higher than that of nitrogen 
replete cultures (initial nitrate concentrations 30 |J,M and 550 
respectively; Figure 2G,H). On day 3, when samples were 
collected for proteome and transcriptional analysis, there was a 
small, but significant (t-test, P<0.05) decrease in the Fv/Fm of 
nitrogen starved cultures compared to nitrogen replete cultures 
(Figure 21). In addition to providing a useful tool for our study, 
these increases in cellular DMSP, in a single microalgal species 
under three diverse growth treatments (osmotic stress, oxidative 
stress and nutrient depletion) highlight the multifunctional role of 
DMSP in the ceU. 

APR Activity 

DMSP is a product of sulphur metabolism and the assimilation 
of this element may therefore limit DMSP synthesis. The 
reduction of 5'-adenylylsulphate (APS), catalysed by the enzyme 
APR, is the key point in regulation of the pathway of sulphate 
assimilation in higher plants [29] and is regulated in a demand 
driven manner by the thiols glutathione and cysteine [30,31], 
amino compounds [32,33], carbohydrates [34,35] and hormones 
[36,37]. To test whether APR might also have a regulatory role in 
the biosynthesis of DMSP in T. pseudonana, its activity was 
measured under increased salinity, increased light intensity and 
nitrogen starvation. 

The APR activity measured in T. pseudonana cultures adjusted 
from a salinity of 10 psu to 35 psu increased slightly, but 
significandy after 48 h (Figure 3A). For the remainder of the 
experiment APR activity was only a litde higher than in the 
cultures kept at 10 psu, despite the substantial increase in 
intracellular DMSP concentration of cultures adjusted to a salinity 
of 35, increasing from below the level of detection to 20 mM 
between day 3 and day 6 (Figure 2B). There was no significant 



difference in cellular protein content between the two salinity 
treatments (Figure S2) There was no diflference in APR activity in 
T. pseudonana cultures exposed to a high light intensity of 
1000 |J.mol m— 2 s— 1 compared to those kept at 50 |j.mol 
m— 2 s— 1 (Figure 3B), whilst again there was a clear increase in 
intracellular DMSP concentration (Figure 2E). Again, no signif- 
icant difference in cellular protein content was observed (Figure 
S2). Under nitrogen starvation APR activity was 2.8-fold higher 
than in nitrogen replete cultures on day 2 (Figure 3C), the first 
time point measured, when an increase in intracellular DMSP had 
already been detected (Figure 2H). However, whilst the DMSP 
concentration of the nitrogen starved cells continued to increase 
throughout the experiment the APR activity gradually decreased 
to the level of the nitrogen replete cultures. As was found in our 
previous study [18] cellular protein content was decreased under 
nitrogen starvation, compared to nitrogen replete cells. The APR 
activities measured in the difierent cultures under normal nitrogen 
supply were between 66 and 209 nmol min— 1 mg— 1 protein 
depending on growth stage and condition of the cultures 
(Figure 3C). This is comparable to the APR activity reported for 
the diatoms Thalassiosira weissflogii and Thalassiosira oceanica which 
ranged from 80 to 200 nmol min— 1 mg— 1 protein [38] and, 
interestingly, two orders of magnitude higher than the activities 
typically found in plants. 

Gao et al. [38] also found that APR activity decreased through 
the growth curve despite a progressive increase in DMSP content 
per cell under nitrogen limitation in the chlorophyte microalga T. 
subcordiformis. In addition, Bochenek et al. [22] showed similar 
decrease in APR activity during the course of batch cultures of the 
haptophyte Emiliania huxkyi, another high DMSP producer. These 
results suggest that there is no clear relationship between 
intracellular DMSP concentration and APR activity in T. 
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Figure 2. Regulation of DMSP concentration. The effect of increased salinity (A,B,C; solid line indicates cultures maintained at 10 psu and 
dashed line indicates cultures transferred to 35 psu), increased light intensity (D,E,F,; solid line indicates cultures maintained at 50 (jmols m^^ s^^ and 
dashed line indicates cultures transferred to 1000 irmols m^^ s^') and nitrogen starvation (G,H,I; solid line indicates nitrogen replete cultures and 
dashed line indicates nitrogen starved cultures) on cell number (A,D,G), the DMSP concentration (B,E,H) and Fv/Fm (C,F,I) of Thalassiosira pseudonana. 
Results are shown as means ± standard deviation from 3 independent cultures. The vertical line indicates the point at which the salinity/light 
intensity of the cultures was adjusted. Asterisks show t-test, P<0.05. Arrows indicate the point at which cultures were sampled for proteomics and 
gene expression analysis. 
doi:1 0.1 371/journal.pone.0094795.g002 
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Figure 3. Regulation of APR activity. The effect of Increased salinity 
(A; solid line, cultures maintained at 10 psu; dashed line, cultures 
transferred to 35 psu), Increased light Intensity (B; solid line, cultures 
maintained at 50 i^mol m"^ s"^ dashed line, cultures transferred to 
1000 |imol m^^ s"') and nitrogen starvation (C; solid line, nitrogen 
replete cultures; dashed line, nitrogen starved cultures) on APR activity 
in Thalassiosira pseudonana. Results are shown as means ± standard 
deviation from 3 independent cultures. The vertical line Indicates the 
point at which the salinity/light Intensity of the cultures was adjusted. 
Asterisks mark significantly different values (P<0.05, T-test). 
dol:1 0.1 371/journal.pone.0094795.g003 



pseudonana. The initial higher APR activity under nitrogen 
starvation might be required for a short term increase in sulphur 
assimilation to initially increase DMSP synthesis, l)ut subserjuendy 
the sulphur for DMSP synthesis might be made available by other 
cell processes. In accordance, the lack of regulation of APR activity 
by increased light intensity suggests that APR is not critical for 
increased DMSP biosynthesis. Given the very high APR activity 
levels it is possible that this is adequate to provide reduced sulphur 
for DMSP synthesis without a need for further up-regulation. Gao 
et al. [38] reported that APR activity of the dinoflagellate 
Heterocapsa triquetra, was only 0.5 to 5 nmol min— 1 mg— 1, which 
is considerably lower than that found for diatoms in that study and 
confirmed here. Thus, as H. triquetra has a high intracellular DMSP 
concentration of approximately 300 mM [39], compared to a 
maximum of 1 5 to 20 mM for T. pseudonana, regulation of APR 
activity may be uncoupled from control of DMSP synthesis. 



Alternatively, there might be a fundamental divergence in the 
regulation of sulphur metabolism between these two groups of 
phytoplankton. This is plausible, as there is evidendy a substantial 
difference in regulation of APR in T. pseudonana and plants. While 
in the diatom APR was not significantly affected by increased 
salinity and was higher in nitrogen starved cells than in controls, in 
Arabidopsis the enzyme is induced by salt and repressed by 
nitrogen starvation [33,40]. 

Expression of Genes involved in Sulphur Assimilation 

Since APR does not seem to be important for control of DMSP 
synthesis, the contribution of other components of sulphur 
metabolism was assessed by comparing transcript levels of genes 
associated with the assimilation of sulphate to cysteine under 
increased salinity, increased light intensity and nitrogen starvation. 
Surprisingly, little coordination was found in the transcript 
responses to the three treatments, including the two isoforms of 
APR (Figure 4). In contrast, in plants these genes are often 
upregulated coordinately e.g. following jasmonate and salt 
treatment, or with the exception of ATP sulphurylase by sulphur 
star\'ation [29]. Transcript levels of APRl (gene ID 3,5690) 
increased under nitrogen starvation, but not with increased 
salinity or light intensity, whereas APR2 (ID 24887) transcript 
levels did increase significantiy with increased salinity, but not with 
the other two treatments. This suggests that different APR 
isoforms may be responsible for the initial increase in APR 
activity seen with nitrogen star\ ati()n and with the increase with 
increased salinity. The lack of changx' in APR transcripts under 
increased light intensity corresponds with the clear lack of 
regulation of APR activity under this treatment. 

Interestingly, the transcript levels of all three serine acetyltrans- 
ferase (SAT, gene IDs 38294, 37497, 16842) genes increased 5-, 
5.7-, and 3.6-fold, respectively, in T. pseudonana under nitrogen 
starvation (Figure 4). In plants SAT is often associated with the 
regulation of sulphur assimilation through the production of O- 
acetylserine (OAS), as modulation of SAT expression affects the 
concentration of cysteine and glutathione [41,42]. Thus OAS 
synthesis might be limiting for cysteine synthesis under nitrogen 
starvation, and consequendy DMSP production in T. pseudonana. 
However, only SATl (ID 38294) expression increased with 
increased salinity and no changes in any SAT isoforms were seen 
with increased light intensity. Therefore, a universal role for SAT 
in control of sulphate assimilation in this diatom species is unlikely. 

The only gene up-regulated under all three growth conditions 
was sulphite rc-duetase (SiR; gene ID 31984; Figure 4), which 
could therefore represent an important point of regulation in 
diatom sulphur assimilation. APR activity in diatoms is about two 
orders of magnitude higher than in plants and could therefore be 
too high to effectively control the flux through sulphate 
assimilation. Thus the control might move to another component 
of the pathway possibly by limiting cysteine synthesis. SiR is a 
good candidate for such control point as reduction of its expression 
in Arabidopsis limits growth [43] . Alternatively, in addition to the 
different responses of APR activity, the general lack of coordina- 
tion in this pathway, across the treatments indicates that sulphur 
assimilation might not limit DMSP production. This is supported 
by our finding that reducing sulphate availability in ESAW 
medium from 25 mM to 5 mM had no effect on DMSP 
concentration in T. pseudonana (Figure SI) and, contrary to what 
has been seen for other marine phytoplankton species, no negative 
effect on growth [44] . With Emiliania huxleyi Bochenek et al. [22] 
found that reducing sulphate concentration to 5 mM reduced 
growth rate by 50% and intracellular DMSP concentration by 
60%. It should be noted that E. huxleyi has an intracellular DMSP 
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Figure 4. Regulation of gene expression. The relative fold change in the transcript levels of genes involved in sulphur assimilation in 
Thalassiosira pseudonana exposed to nitrogen starvation, increased salinity and increased light intensity quantified by qRT-PCR. Gene abbreviations 
are as follows: Sulphate transporter (Trans), ATP sulphurylase (ATPS), APS reductase (APR), OAS thiollyase (OASTL), Serine acetyltransferase (SAT), 
Sulphite reductase (SIR). Numbers denote different isoforms of the enzymes (for gene IDs see Table SI). Results are shown as means ± standard 
deviation from 3 independent cultures. Asterisks mark values significantly different between treatments and control at P<0.05 (T-test). 
doi:l 0.1 371 /journal.pone.0094795.g004 



concentration several fold higher than that of T. pseudonana and so 
is likely to have higher sulphur requirements than the diatom [45] . 
To add complexity to the regulation of sulphur assimilation and 
DMSP synthesis, at a number of steps in the pathway there are 
multiple isoforms that can catalyse the reaction, and these might 
well have different kinetic properties and/or be localised to 
different cellular compartments. The difiFerential regulation of 
isoforms seen here thus warrants further investigation. 



concentration. Although the fold change was not greater than 
1.5-fold in all cases, at least one spot in the chain increased by 1.5- 
fold under each treatment. These spots were interesting because 
they could be same protein possibly separated due to posttrans- 
lational modifications that affect the proteins isoelectric point. Five 
spots were picked and indeed all identified as a phosphoenolpyr- 
uvate carboxylase (PEPC; ProdD 268546; table 1). 



Identification of Proteins associated with DMSP synthesis 

In an unbiased approach to identify enzymes involved in DMSP 
biosynthesis, we identified changes in the proteomes of T. 
pseudonana cells in response to increased salinity and light intensity 
and compared these to our previous dataset on the response of T. 
pseudonana to nitrogen starvation [18]. The separation of protein 
extracts by 2-dimensional gel electrophoresis yielded 3310 
distinguishable protein spots. There were 479 spots that changed 
by more than 1.5-fold in relative abundance with q<0.05 
(determined by t-test and corrected by FDR) under any of the 
three treatments. We had hypothesised that any protein found to 
change in abundance in the same direction under all three 
treatments would potentially be involved in DMSP biosynthesis, 
however we found that most of the changes were treatment 
specific, with very few proteins changing in abundance under 
multiple treatments. Only 1 spot changed under all three 
treatments, however, the 1,5-fold cut-off was met only in the 
nitrogen and high light dataset. (Figure 5). This spot (ID: 00569) 
belongs to a horizontal chain of spots that all increased under all 
three growth conditions that increased intracellular DMSP 




Figure 5. Global proteome changes. Venn diagram showing the 
number of protein spots Identified by 2-dlmenslonal gel electrophore- 
sis, changing In relative abundance by more than 1.5-fold (t-test, FDR, 
q<0.05) In Thalassiosira pseudonana under nitrogen starvation. 
Increased salinity and Increased light Intensity. 
doi:1 0.1 371/journal.pone.0094795.g005 
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Table 1. Proteins related to sulphur assimilation, methionine metabolism, and photosynthesis changing in abundance under 
nitroaen starvation (N), increased salinity (S) or increased liaht intensity (L). 





Protein ID 


Protein Name 


Fold Change 




FDR (q) 










N 


S 


L 


N 


S 


L 


260934 


BrsnchEd-chsin-smino-scid sminotrsnsfcrssG 


6.4 


1.3 


— 1.2 


0.0062 


0.2804 


0.2638 


270365 


Sulphite reductase'^'''^ 


1.6 


1.2 


1 .1 


0.01 87 


0.2804 


0.2998 


20797 


^-3Hpnnc\/mP>triirmin,i f'sAh^^- Hpin^mHonl" 
J aut^i y 1 1 1 1 C LI 1 1 U 1 1 1 1 it^ LJmvi/ Ut^LJt^l lU til 1 L 

Methyltransferase"'^' 


-1.7 


3.6 


1 2 


0 0516 


0 0559 


0 1854 






2 0 


3 6 


1.3 


0.0492 


0.0574 


0.1156 


644 


Adenosine kinase 


1 1 


2 5 


1.2 


0.2455 


0.0179 


0.0558 


21815 


Adenosylmethionine synthetase 


14 


2 9 


1.1 


0.0469 


0.0277 


0.25564 


27273 


Methylenetetrahydrofolate reductase 


1 Q 


2 0 


1.5 


0.4515 


0.0582 


0.0257 


28496 


Adenosylhomocysteinase 


1 y 


2 6 


1.3 


0.021 


0.0569 


0.0765 


25402 


Fucoxanthin chl a/c iight-liarvesting protein 


_ 1 Q 


1 2 


-3.6 


0.424 


0.3975 


0.0118 


22747 


Fucoxanthin chlorophyll a/c light-harvesting protein 


-1.0 


1.0 


-3.5 


0.4793 


0.4825 


0.0127 


bd1048 


Photosystem 1 iron-sulphur centre*^^ 


1.0 


-1.2 


-2.6 


0.3337 


0.0957 


0.0159 


30385 


Fucoxanthin chlorophyll a/c protein-LI818 clad 


-1.2 


-1.2 


-2.5 


0.2049 


0.3020 


0.0127 


24080 


Fucoxanthin chl a/c light-harvesting protein 


-1.1 


1.2 


-2.4 


0.346 


0.405 


0.0285 


38583 


Fucoxanthin chlorophyll a/c protein 


1.0 


-1.4 


-2.4 


0.4518 


0.0408 


0.0118 


268546 


Phosphoenolpyruvate carboxylase 


2.5 


1.7 


1.4 


0.0161 


0.1789 


0.0547 






2.0 


1.5 


1.4 


0.0499 


0.2037 


0.0473 






1.8 


1.3 


1.4 


0.0208 


0.1746 


0.0244 






1.6" 


1.7" 


1.5" 


0.023 


0.0661 


0.0801 






2.5 


1.2 


1.6 


0.0062 


0.0408 


0.0181 



"^Manual Annotation. 

''Supported by BlastP {E<3x10""). 

*^Supported by conserved domains identified through Pfam. 
"Combined with another protein making fold change imprecise. 

Protein names are based on UniProtKB unless otherwise stated and Protein IDs are from the Joint Genome Institute T. pseudonana genome version 3 (http://genome.jgi- 
psforg/Thaps3/Thaps3. home. html}. 
doi:l 0.1 371/journal.pone.0094795.t001 



Since the response of the T. pseudonana proteome was quite 
different between the different treatments we picked the spots with 
the strongest changes from each treatment. For the sahnity and 
hght treatments the 10 spots with the greatest fold change, and 
sufficient quantity were picked from the gel, and in the nitrogen 
starvation study aU spots that changed were picked [18]. In total 
135 spots were picked, from which 16 could not be identified, 
resulting in identification of 84 unique proteins, by MALDI-TOF 
MS analysis of the tryptic digests, that changed in relative 
abundance under one or more treatment (Table S2), including 
proteins related to sulphur metabolism. It is important to note that 
the absence of a protein among the analysed spots does not prove 
its stable abundance as not all regulated proteins can be identified 
by 2-dimensional electrophoresis. For example gels are not 
optimised for membrane bound proteins and very low abundance 
proteins might also not be detected. 

Increased Salinity. Among proteins which increased in 
abundance specifically with increased salinity (Tables 1, S2), the 
proteins involved in the active methyl cycle (Figure 6), which 
salvages methionine from methyl transferase reactions, were most 
prominent. An S-adenosylmethionine (SAM) synthetase (ProtID 
21815), which catalyses the formation of SAM from methionine, 
increased by 2.9-fold. SAM is used as a methyl donor in a wide 
variety of methyltransferase reactions and accordingly a 3.6-fold 
increase in the abundance of a SAM-dependent methyltransferase 
(ProtID 20797) was measured. Enzymes of this family have a 



diverse range of functions [46] including synthesis of various 
metabolites. The SAM-dependent methyltransferase identified 
(ProtID 20797) has similarities to a sarcosine/dimethylglycine 
methyltransferase (BlastP E = 3x10 — 32) which provides an 
alternative route for synthesis of the osmolyte glycine betaine 
(GBT) through the methylation of glycine. This metabolic route of 
GBT synthesis is distinct from the pathway of choline oxidation 
found in higher plants and it has been identified in a number of 
halotolerant bacteria and cyanobacteria [47,48] and the red alga 
Galdieria sulphurmia [49]. 

The S-adenosylhomocysteine produced by such methyltrans- 
ferase activity can be recycled to homocysteine through the action 
of adenosylhomocysteinase; this enzyme (ProtID 28496) increased 
2.6-fold with increased salinity. Adenosine is the second product of 
this reaction and, correspondingly, an adenosine kinase (ProtID 
644), which catalyses the reversible phosphorylation of adenosine 
to adenosine monophosphate, increased 2.5-fold. This enzyme 
might remove adenosine, thereby promoting the activity of 
homocysteinase. Together the enzymes of the active methyl cycle 
have the potential to influence the free methionine pool and 
thereby its availability for DMSP synthesis under increased 
salinity. An increase in the abundance of enzymes associated with 
the active methyl cycle with increased salinity has also been 
reported in the sea-ice diatom Fragilariopsis cylindms [50]. 

The abundance of methylenetetrahydrofolate reductase (ProtID 
27273) increased 2-fold under increased salinity. This enzyme 
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Figure 6. A scheme of the active methyl cycle. Enzymes that increased in abundance with increased salinity In T. pseudonana are marked with 
bold arrows. 
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reduces 5,10-methylenetetrafolate to 5-methyltetrahydropteroyl 
tri-L-glutamate, which donates its methyl group for synthesis of 
methionine from homocysteine. The disruption of this enzyme in 
the bacterium Streptomyces lividans leads to methionine auxotrophy 
[51]. Methylenetetrahydrofolate reducta.se might therefore be 
limiting methionine .synthesis and DMSP production in T. 
pseudonana. 

Increased Light Intensity. Methylenetetrahydrofolate re- 
ductase (ProtID 27273) increased 1.5-fold also under increased 
light inten.sity, suggesting that there may be some similarity in the 
response of sulphur metabolism to increased salinity and increased 
hght intensity. However, enzymes of the active methyl cycle 
remained unchanged. 

The most obvious response of T. pseudonana to increased light 
intensity is the degradation of photosynthetic proteins; a weU- 
delined response in photosynthetic organisms [52-55]. Five light 
harvesting proteins (ProtlDs 25402, 22747, 30385, 24080, 38583) 
and one photosystem I centre (ProtID bdl563) decreased in 
abundance relative to the control (Table 1), suggesting that this is 
also the case in the response of T. pseudonana to increased light 
intensity, although reduced translation or transcription cannot be 
ruled out. The degradation of photosynthetic proteins could 
increase the availability of amino acids such as cysteine and 
methionine for DMSP synthesis. Grone and Kirst [56] first 
demonstrated that protease inhibitor delayed DMSP accumulation 
in the Prasinophyte Tetraselmis (Platymonas) subcordiformis suggesting 
that protein degradation has a role in the regulation of DMSP 
synthesis. Futhermore, Lyonn et al [50] have also proposed that 
degradation of light harvesting proteins could be a potential source 
of cysteine and methionine for DMSP induction based on their 
study of the F. cylindrus proteome response to salinity. 

Nitrogen Starvation. Under nitrogen starvation a decrease 
in cellular protein level was measured in T. pseudonana and there 
was also evidence of increased protein catabolism [18], which 
would influence the amino acid pools in the same way as the 
degradation of photosynthetic proteins under increased light 
intensity. On the other hand, the proteome response to nitrogen 
starvation was quite different to that of increased salinity. The 
abundance of methylenetetrahydrofolate reductase was un- 
changed suggesting that this enzyme does not contribute to the 
increase in DMSP measured under this treatment. In addition the 



abundance of many enzymes of the active methyl cycle were 
reduced under nitrogen starvation, with the SAM-dependent 
methyltransferase identified in this study decreased by 2.1 -fold. 
This is not entirely surprising since a reciprocal relationship that is 
dependent on nitrogen availability has been found between 
cellular DMSP and the nitrogen based osmolyte, glycine betaine 
(GBT) in continuous cultures of T. pseudonana [57]. However, this 
again suggests that, whilst the active methyl cycle might have a 
role in increasing the free methionine pool and DMSP synthesis 
with increased salinity, this is not the case regarding nitrogen 
starvation. 

The protein with the greatest increase in abundance identified 
in this study was a branched chain aminotransferase (ProtID 
260934) that increased 6.4-fold under nitrogen starvation, 
however it did not change under increased salinity or light 
intensity. This enzyme might be a candidate for the first step of 
DMSP biosynthesis, catalysing the production of MTOB from 
methionine by transamination. However, further research would 
be required to confirm this, such as substrate specificity testing and 
demonstrating that over expression or silencing of the gene has an 
effect on DMSP synthesis. If this enzyme was found to catalyse the 
transamination of Methionine to MTOB it would raise questions 
as to why it did not increase with increased salinity or light 
intensity. 

The only enzyme found to increase in abundance by more than 
1.5-fold under all three treatments that increased the intracellular 
DMSP concentration of T. pseudonana was a PEPC (ProtID 
268546). Oxaloacetate, produced by PEPC, is used in the synthesis 
of aspartate, which is a precursor of methionine. O-phosphoho- 
moserine (OPH; Figure 1) is an intermediate between aspartate 
and methionine, and is the substrate of both cystathionine a- 
synthase (CgS), forming cystathionine, and threonine synthase, 
forming threonine. The balance of these processes is an important 
point of regulation in methionine synthesis in higher plants [58]. It 
is plausible that increased PEPC activity is a response to increased 
demand for aspartate and DMSP synthesis. In plants, particularly 
those with C4 photosynthesis, PEPC is also induced by salt stress, 
it is linked to pH homeostasis and carbon/nitrogen balance [59]. 
Similar roles in diatoms cannot be excluded, indeed, the increased 
expression of PEPC and aminotransferases can be used for 
intracellular nitrogen recycling to maintain C/N balance [60]. As 
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T. pseudomm expresses all enzymes of the C4 photosynthesis cycle 
[61] it can be speculated that the increase in PEPC might serve as 
a CO2 pump to reduce photorespiration and loss of C and N. 
Should these be the mechanisms of action of PEPC, the increase 
DMSP might be just a mechanism of dissipation of excess sulphur. 
However, this would most probably trigger a repression of 
sulphate uptake and assimilation, which was not the case. 

What does this tell us about the regulation of DMSP 
production? 

The lack of coordination in transcript level throughout the 
pathway of sulphur assimilation in these experiments suggests that 
sulphur does not limit DMSP production in T. pseudonam, as is the 
case for sulphur compounds in higher plants. This would not be 
entirely surprising since marine algae evolved in a very different 
sulphur environment to terrestrial plants [44]. The ocean is a 
sulphur reservoir where this nutrient is unlikely to limit 
productivity, whereas in terrestrial habitats sulphur availability 
does limit growth. The APR activities measured here in T. 
pseudonana, and in other algal species [38], are more than two 
orders of magnitude higher than in plants. Gao et al. [38] suggest 
that these algal species have a high capacity to assimilate sulphate 
and may build up intracellular stores that are used in DMSP 
production rather than newly assimilated sulphate. 

Another possible reason for the lack of coordination in the 
proteome response of T.p.seudonana under conditions that increase 
in DMSP synthesis is that this process might not depend on the 
synthesis of new proteins. There are many possible levels of 
regulation, so that possibly aUosteric regulation(s) may control 
DMSP synthesis, as is the case for methionine synthesis. 
Threonine synthase, which competes with CgS for OPH, is under 
allosteric control and this regulates the flow of carbon into 
threonine' or methionine. 

There are some patterns in the transcript and protein responses 
to the treatments used here that could highlight potential points of 
regulation other than sulphate assimilation. Interestingly, a 
number of these changes occurred in the branches that bring 
carbon and nitrogen skeletons to the central pathway of sulphur 
assimilation. These include increased transcript levels of SAT 
under nitrogen starvation, increased methyltetrahydrofolate re- 
ductase protein abundance under increased saJinity and light 
intensity, and increase PEPC protein abundance under all three 
growth conditions. This may suggest that carbon and/ or nitrogen 
skeletons limit DMSP synthesis rather than sulphur. Only 
transcript levels of SiR were increased under all three growth 
conditions suggesting that this enzyme might be a key controlling 
point of sulphur assimilation in diatoms rather than APR as in 
higher plants. 

In a study that focussed on the effect of increased salinity on 
DMSP synthesis in the sea-ice diatom Fragilariopsis cylindms, Lyon 
et al. [50] identified candidate enzymes for the four classes of 
enzymes that Gage et al. [14] proposed to catalyse the conversion 
of methionine to DMSP (aminotransferase, reductase, methyl 
transferase, decarboxylase). Although their activity has yet to be 
confirmed, it is interesting that using a comparable proteomic 
approach, but looking f )r proteins that changed in three different 
tr(;atmc'nts that increase inrra( c-llular DMSP concentration, these 
proteins were not identified. This could be due to differential 
regulation between these species, indeed these diatoms inhabit 
quite different environments. It is also possible that these proteins 
did increase under increased salinity, but were not among the 
highest changing and therefore not picked, since not all protein 
spots were analysed in this study. However, neither were they 
identified in our nitrogen starvation treatment where all the 



proteins that changed significantly were analysed. It is possible that 
the changes in protein abundance observed by [50] are specific to 
increased DMSP content with increased salinity and that 

alternative processes might be responsible for the increase in 
DMSP under nitrogen starvation. There was a very high 
consistency between the 2D gels analysed in our study, however, 
we cannot dismiss the possibility that an increased number of 
replicates might lead to identification of additional significant 
differences. 

It is thought that DMSP has multiple functions within the ceU, 
and indeed, in the experiments here it might be considered to be 
used as an osmolyte, antioxidant and possibly as an overflow 
metabolite. Perhaps the lack of coordination in the response of the 
pathway of sulphur assimilation to these growth conditions 
indicates that different regulatory mechanisms, each potentially 
leading to an increase in the methionine pool, might relate to its 
different roles in the ceU under the different conditions. Under 
increased salinity the production of glycine betaine via the 
methylation of glycine, might affect the methionine pool through 
the active methyl cycle; of which many of the enzymes were seen 
to increase in abundance. The increase in the abundance of 
methylenetetrahydrofolate reductase under increased salinity and 
light intensity is likely to increase 5-methyltetrahydropteroyl tri-L- 
glutamate supply for methionine biosynthesis, thereby increasing 
the availability of this amino acid. Methionine availability could be 
increased under increased light intensity and nitrogen starvation 
by protein degradation. Perhaps there is no individual limiting step 
or 'on-switch' for DMSP synthesis, instead under different 
conditions different components could be limiting. This might 
depend on the carbon and nitrogen status of the cell, each of 
which win be different under the three conditions tested here. 
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